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An embedded selective recoupling method is proposed which is base d on the ide a of embedding the 
recently proposed deterministic selective recoupling scheme of lYamaguchi et all Q] into a stochas- 
tic dynamical decoupling method, such as the recently proposed Pauli-random-error-correction- 
(PAREC) scheme 2j. The recoupling scheme enables the implementation of elementary quantum 
gates in a quantum information processor by partial suppression of the unwanted interactions. The 
random dynamical decoupling method cancels a significant part of the residual interactions. Thus 
the time scale of reliable quantum computation is increased significantly. Numerical simulations are 
presented for a conditional two-qubit swap gate and for a complex iterative quantum algorithm. 
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I. INTRODUCTION 

Stabilizing quantum systems against uncontrolled per- 
turbations is one of the major challenges in quantum 
information science. In comparison with measurement- 
based quantum error correction, at first sight dynamical 
decoupling methods seem to suffer from the disadvan- 
tage of not being capable of correcting unwanted per- 
turbations perfectly. However, typically measurement- 
based quantum error correction requires ancillary quan- 
tum systems and complicated measurement and recovery 
operations |3| which are difficult to realize by nowadays 
technology. By contrast dynamical decoupling methods 
are able to suppress unwanted perturbations significantly 
even over long interaction times without requiring ad- 
ditional ancillary quantum systems Thus, 
in view of nowadays experimental difficulties in control- 
ling many-particle quantum systems these latter methods 
offer interesting perspectives, in particular, for the first 
generation of quantum information processors. 

Recently a deterministic dynamical recoupling scheme 
was proposed for dipole-coupled nuclear spins in a crys- 
talline solid m. While in all previously proposed sim- 
ilar schemes 0, IS U the evolution-time overhead 
grows linearly with the number of spins, this particu- 
lar scheme leads to an evolution-time overhead which 
is independent of the number of spins involved. Thus, 
it appears to be well suited for the stabilization of 
quantum information processors against unwanted inter- 
qubit interactions. This recoupling scheme uses par- 
ticular combinations of fast broadband and slower se- 
lective radio-frequency fields to turn off all couplings 
except those between two particularly selected ensem- 
bles of spins. Thereby, spins within each ensemble rep- 
resenting a particular logical qubit are decoupled [Tlj. 
Furthermore, cross-couplings between selected ensembles 
are avoided by requiring that qubit couplings have to 
be much stronger than any other couplings within each 
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ensemble. Unwanted couplings are suppressed up to 
second-order average Hamiltonian theory with the help 
of time-symmetric pulse sequences. Despite many ad- 
vantages in this recoupling scheme the residual higher- 
order interactions accumulate coherently thus leading to 
a quadratic-in-time decay of the fidelity of any quantum 
state 12|. This restricts the achievable time scales of 
reliable quantum computation significantly. 



In this paper it is demonstrated that the performance 
of this recoupling scheme can be improved significantly 
by embedding it into a stochastic decoupling scheme. For 
this purpose the recently proposed Pauli-random-error- 
correction (PAREC) method Q turns out to be partic- 
ularly useful. Thereby, the stochastic decoupling scheme 
destroys the coherent accumulation of higher-order resid- 
ual interactions to a large extent so that the fidelity decay 
of any quantum state is slowed down significantly to an 
almost linear-in-time one. As a result reliable quantum 
computation can be performed on significantly longer 
time scales. 



This paper is organized as follows: The basic ideas un- 
derlying the recently proposed deterministic recoupling 
scheme of Ref. Q are summarized briefly in Sec. Ill Al for 
the sake of completeness. In Sec. Ill Bl a simple restricted 
embedded decoupling scheme is introduced. Though it 
already leads to first improvements in comparison with 
the deterministic selective recoupling scheme of Ref. |l[ , 
its error suppressing properties can still be improved sig- 
nificantly by an additional simple symmetrization pro- 
cedure. In Sec. IIIII the stabilization properties of this 
symmetrized embedded recoupling scheme are analyzed 
for a unitary conditional two-qubit swap gate. In Sec. 
IIVI its stabilizing properties are investigated by apply- 
ing it to the iterated quantum algorithm of the quantum 
sawtooth map . 
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II. SELECTIVE RECOUPLING BY EMBEDDED 
DYNAMICAL DECOUPLING 

A. Deterministic selective recoupling of qubits 

In this section the basic ideas underlying the recently 
proposed recoupling scheme of Ref. are summarized. 
In particular, the form and magnitude of the residual 
higher-order interaction is discussed which cannot be 
suppressed by the suggested pulse sequences. 

Let us consider n q nuclear spin- 1/2 systems in a crys- 
talline solid which are interacting with an external static 
magnetic field in z-dircction. In the rotating wave ap- 
proximation their Hamiltonian is given by 

n q -l n„-l 

k=0 k<l=0 
= Hz + Hd 

(1) 

with the Pauli spin operators a x , a y , and a z . Thereby, 
the Larmor frequencies uj^ of the first term characterize 
the interaction strengths of these spins with an external 
inhomogeneous magnetic field so that these spins can be 
addressed individually. The second term of the Hamilto- 
nian |JQ describes the dipole-dipole interaction of the nu- 
clear spins with the coupling strength J^i between spins 
k and I being inversely proportional to the cubic power 
of their distance. 

If these nuclear spins are used as qubits of a quantum 
memory, for example, one has to protect them against the 
perturbing influence of the interaction Hamiltonian 10 . 
In the framework of a deterministic decoupling scheme |6( 
this may be achieved by an appropriate sequence of fast 
electromagnetic pulses. Thereby, a series of fast global 
7r-pulses is applied in order to suppress the Zeeman term 
Hz while leaving the dipole-dipole coupling term Hd in- 
variant. This latter term can be suppressed by the well 
known WHH-4 pulse sequence This sequence con- 

sists of four fast it/ 2- pulses applied at times r, 2r, 4r 
and 5t. Thus, the resulting unitary time evolution after 
this pulse sequence, i. e. at time T c = 6r, is given by 

6 

U(t = 6r) = L> 1 t (-r)LT 1 (r) = TT e -**W» 

f=i (2) 

U x {t) = e- lilT / h Pye- lflT/h P x e- lflT/h 

with the interaction-picture (toggled) Hamiltonians 
Hi = H 6 = H D , H 2 = H 5 = P X H D P X and H 3 = 
H4 = P x P y HDPyP x - Here, the unitary transforma- 
tion P x results from a fast global 7r/2-pulse witch gen- 
erates rotations around the x-direction, i. e. P x = 
(^^Iq 1 exp(— ia x n/4) = PI. As a consequence, in 
zeroth-order average Hamiltonian theory (AHT) 6] the 



time-averaged Hamiltonian vanishes, i. e. 

_s_(0) 1 JL - 

Due to the time reversal symmetry of the WHH-4 pulse 
sequence, i. e. H(t) = H(T C -t), in AHT all odd higher- 

_^(2n+l) 

order Hamiltonians vanish, i. e. H D = for n € No- 
li these nuclear spins are used as qubits of a quantum 
information processor one also has to implement univer- 
sal sets of unitary quantum gates. In particular, one 
needs to be able to implement two-qubit entanglement 
gates, such as controlled phase gates. This can be accom- 
plished by recoupling qubits selectively with the help of a 
Super- WHH pulse sequence as proposed in Ref. . Such 
a Super- WHH sequence recoupling qubits k and I consists 
of three WHH-4 sequences applied to the toggled Hamil- 
tonians Hf z = a k a l z H D a k d l zl Hf y = a k & l y H D & k & y , and 

H kl x — a k a l x H D a k a l x , respectively. Correspondingly, 
there are 18 time periods of duration r during which the 
time evolution is described by the double-toggled Hamil- 
tonians Hi = H kl z , H 2 = P x H kl z P x , etc.. The appro- 
priate WHH-4 pulse sequence of the H kl y Hamiltonian, 
for example, is illustrated in Fig. ^ where free evolution 
denotes the time evolution according to the Hamiltonian 
Hd over a time interval of duration r. The quantum 
gates resulting from these WHH-4 sequences are denoted 
by W xy , W kl z , and W k x , respectively . The Super- WHH 
sequence is finally obtained from a combination of these 
latter quantum gates preceeded by the corresponding 
time reversed sequence (compare with the inner part of 
Fig. |2J. As a consequence pj, this Super- WHH sequence 
yields the zeroth-order time-averaged Hamiltonian 

~ (0) 1 ^ ~ 

^ = ^ E % = J ^ + *A + W w 
3=1 

with the renormalized zeroth-order recoupling strength 
jjS — (Jjy/4) x (8/9). Due to the time reversal symme- 
try of the Super- WHH sequence, in AHT all odd-valued 
higher order recoupled Hamiltonians vanish, i. c. i. c. 

H D = for n E N . 

With the help of the zeroth-order recoupled Hamilto- 

— (0) rn 

nian H D of Eq. J3J one can approximate unitary two- 
qubit quantum gates of the form 

llf = exp(-i(&*dl + & k v o- l y + a k a l z )4>). (5) 

Thereby, for a particular value of the phase <f> one has 
to adjust the time r between two successive pulses of a 
WHH-4 sequence and the number of times n sw hh a Super- 
WHH sequence has to be applied according to the rela- 
tion 

J$ • n swhh 36T/h = 4> (6) 



FIG. 1: Schematic representation of the unitary W*y quantum gate acting on qubits k and I: Free evolution indicates time 
evolution according to the Hamiltonian Ho over a time interval of duration r. 
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FIG. 2: Schematic representation of a conditional two-qubit gate (t/^-gate) obtained by recoupling qubits k and I according to 
Eq. The W 7 * gates are obtained by reversing the order of the broadband pulses suppressing the Zeeman term. Residual 
second-order terms of AHT can be eliminated by the restricted randomization step accomplished by random selective 7r-pulses 
ov . . . Thereby n t k has to be equal to Vi : i to ensure that the wanted gate action is not disturbed. Still remaining terms are 
symmetrized by random 7r/2-pulses R ai (n/2) — exp(— i<r ai 7r/4) , on € {x,y,z}. Either condition (|SJ or condition 11 H has to 
be fulfilled depending on whether the original Super- WHH or the symmetrized Super- WHH sequence is used. 



(compare with Fig. [5J. However, because of the residual by the Super- WHH pulse sequence this implementation 
higher-order interactions which have not been canceled of a two-qubit quantum gate gate is only approximate. 

I 

The error resulting from these residual higher-order interactions is dominated by the second-order term of AHT 
which is generally given by 

^( 2 ) 1 f Tc f* 3 f t2 / - ~ - \ 

H =-vfJ o dt 3j dt 2j o *i([ J ff(^),[ff(i 2 ),i/(ii)]] + [[ J ff(t3), J ff(i2)],ff(ti)])/ft 2 (7) 

for an arbitrary time-dependent Hamiltonian H(t). In the case of the Super- WHH sequence we have H{t) = Hj if 
(j — 1)t < t < jr and T c — 36r. Therefore, according to AHT the lowest-order correction to the recoupled Hamiltonian 
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of Eq. J2J is given by 

t D = [ & * a * (~ 322 J aVafe + 446 J a \.J a , + 3628 J a/ J afe J fci - 2906 - 1370 J 2 ,J fci 



+308J a V afe + 308 J a \J a/ + 3208 J a iJ ak J kl - 2588 J 2 fe J fci - 2588 J 2 ^ 



(8) 

o*St (+446 J a y afe - 322 J a \. J a/ + 3580 J al J ak J k i - 1922 J 2 fc J H - 345 8 J 2 ^)] T 2 /(fr 2 ■ 1728) + . . . . 



Thereby, only terms of the form a\a 1 ^ = \®---®\®b\®\®---®\®o l p®\®---®\, (a,j3 ^ {x,y, z}) are indicated 
as all other terms are irrelevant for our subsequent discussion. 



As a consequence, the Hamiltonian resulting from recou- 
pling of qubits k and I by a Super- WHH sequence is of 
the form 

H%' = 4? (Z k A + o\a y + a k z a l z ) + O (J(Jr/hf) . (9) 

To estimate the resulting error affecting the unitary 
gate U¥ generated by fljj we have to multiply the 
O (J(Jr/h) 2 ) term of © by the duration n swhh 36r of 
this quantum gate. Using condition © we see that for a 
fixed value of the phase 4> the resulting error depends on 
™ S whh according to J 3 r 2 • n swhh 36r/?i ~ 3 /(36n swhh ) 2 . 
Therefore, the resulting contribution to the fidelity fu* 
of any quantum state \^} is given by 

fu> = muiuLm = 1 - O(0 6 /(36n sw hh) 4 ). (10) 



B. Randomization and selective recoupling 

The residual interaction described by the Hamiltonian 
iJSJ can be suppressed significantly by embedding the re- 
coupling scheme of Sec. Ill Al into a stochastic decoupling 
scheme, such as the recently proposed PAREC scheme 
. For this purpose we choose at random an n 9 -fold ten- 
sor product of Pauli-matrices a ri <E> <J Ti x ® • • • ® oy< _ x 
(r G {0, x, y, z}, (To = 1) and apply it before and after the 
i-th Super- WHH sequence. This way each deterministic 
Super- WHH sequence is embedded within two statisti- 
cally independent random Pauli operations. In contrast 
to a usual dynamical decoupling scenario in our case 
we have to choose the Pauli-matrices in such a way that 

they leave the ideally recoupled gate Hamiltonian H D 
of Eq. (0J invariant. This can be achieved by impos- 
ing the restriction that the randomly chosen statistically 
independent Pauli spin operators have to be identical 
for qubits k and I for each Super- WHH sequence, i. e. 
5y 4 k = & ri , . This restriction assures that terms of the 

form o\cr a in H^ 1 remain invariant (compare with Fig. 

-£-(2) 

0). Since H D contains no terms of the form ffj^ with 
a 7^ (3 (compare with Eq. (jHJ) the Pauli-matrices for 
qubits k and I can always be omitted, i. e. chosen to be 
the identity, r i:k = r it i = 0. 



— ( 2 ) 

The only terms of the Hamiltonian H D which cannot 
be eliminated by this constrained randomization method 
are the ones containing terms of the form <j\o l a (a £ 
{x,y,z}) which are shown in Eq. (JSJ. However, by an 
additional symmetrization these terms can be made rota- 
tionally invariant so that they can be cast into the form of 
Eq. (@J. Thus, for a given value of (j> these terms lead to a 
renormalization of the values of the required gate param- 
eters r and n sw hh- This rotational symmetrization can 
be achieved by selective 7r/2-pulses which induce unitary 
transformations of the form R a (ip) — exp(—ia a /2 ■ tp) 
with a 6 {x,y, z}. For this purpose one chooses one of 
the three unitary transformations R*. (ir/2)R! a . (it/ 2) act- 
ing on qubits k and / randomly and applies it before and 
the corresponding inverse transformation after the i-th 
Super- WHH sequence (compare with Fig. |2J). This way 
the coefficients of the cr„<7^-terms are permuted in the 
relevant toggled Hamiltonians. As a consequence one ob- 
tains the statistically and rotationally averaged second- 
order contribution 

tl = (a k x a l x + a k a l y + a k a l z ).J^r 2 /h 2 (11) 

with 

(2) ^ ( 1 

^kl ~ \ ~^(JalJak + Ja k Jal) 

a 

+ J^JalJ ak J k l ~ -^-(j£ k Jkl + JllJkl^j- (12) 

By this combined randomization and symmetrization 
method the improved recoupled Hamiltonian 

Hf = {j^ + J^(r/nf + 0{j(Jr/hf)). 

(cr k x cr l x +a k y a y +a k a l z ) (13) 

is obtained. In contrast to Eq. JJjJ, now the effective re- 
coupling strength is renormalized and the residual error 
is suppressed up to fourth-order in the small coupling 
parameter Jt/H <C 1. Thus, in order to implement a 
L^-gate, for example, we now have to choose the renor- 
malized characteristic parameter r' in such a way that 
the condition 

(4? + 4? (t'A) 2 ) ^swhh36r' = (14) 



is fulfilled. As a result, in general the required time of 
free evolution r' depends on the chosen qubit pair (k, I). 

III. ERROR ANALYSIS OF A DYNAMICALLY 
DECOUPLED UNGATE 

In this section the stabilizing properties of selective re- 
coupling by the embedded dynamical decoupling method 
of Sec. IIIBl is investigated for a unitary phase gate as de- 
scribed by Eq. JSJ). As shown in Eq. lfTU|l the fidelity of 
a unitary phase gate which is realized by recoupling 
of qubits k and I with the help of the average Hamilto- 
nian of Eq. @ deviates from unity by terms of the or- 
der of O (</f> 6 /(36n sw hh) 4 ) • Thereby, n sw hh and r denote 
the number of required iterations of the Super- WHH se- 
quence and the time required for the intermediate free 
evolution as determined by relation JSJl, respectively. 

In order to estimate the improvement achievable with 
the help of the embedded recoupling scheme of Sec. Ill Bl 
we start from the observation that in the case of a time 
independent Hamiltonian H a rigorous lower bound of 
the mean fidelity E/(T) of a stochastic dynamical de- 
coupling scheme is given by 0, 0| 

E/(T) > 1 - \\H\\ 2 AtT/h 2 (15) 

for \\H\\ 2 AtT/h 2 <C 1. Thereby, the norm of the Hamil- 
tonian \\H\\ is defined by its largest eigenvalue, At is 
the time interval between successively applied uncorre- 
cted random pulses, E denotes statistical averaging, and 
T is the interaction time. This result can be used to 
estimate the mean error bound of the fidelity of the cor- 
responding embedded decoupling scheme in which the 
recoupling procedure leading to the Hamiltonian © is 
embedded into statistically independent random Pauli 
operations. For this purpose we make the replacements 
0(\\H\\) -> J(Jr/h) 2 , At -> 36t, T n swh h36r. From 
relation JJJJ we obtain the mean error estimate of the 
fidelity to be of the order of O 07(36 4 n|? whh )) . 

In Fig. (bottom) the mean fidelity /swap °f a um ~ 
tary U w /4-gate and its dependence on the number of per- 
formed Super- WHH sequences n sw hh is depicted. In these 
numerical simulations this unitary quantum gate is real- 
ized by recoupling of the two central qubits of a linear 
four-qubit chain. Apart from an irrelevant global phase 
this unitary U n /4-gate is nothing but a Swap-gate (com- 
pare with Fig. @J. The mean fidelity of Fig. [3] was ob- 
tained by averaging over all 2 4 orthonormal initial states 
of the computational basis. The statistical averaging was 
performed over 100 runs with statistically independent 
realizations of the random pulses involved. The coupling 
strength is assumed to be constant for adjacent qubits 
and to be vanishing between all other qubits (compare 
with Fig. [3] (top). Fig. [3] (bottom) demonstrates that 
the fidelity (diamonds) resulting from non-randomized 
Super- WHH pulse sequences can be fitted well by a func- 
tion of the form exp(— c/n 4 hll ) with c « 0.20. This is 
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FIG. 3: The fidelity (bottom) of the *7*/ 4 -gate on a lin- 
ear four-qubit chain (top) as a function of the number 
of repetitions n sw hh of the Super-WHH scheme: original 
Super- WHH sequence (diamonds), unsymmetrized embedded 
scheme (circles), and symmetrized embedded scheme with 
adapted pulse interval r' according to Eq. I|14fl (squares). 
The solid lines represent the fitting functions exp(— c/n* whh ) 
and exp(— c rs /rig whh ) with c = 0.2 and c rs = 0.41. 

consistent with the simple estimate (|10fl . Using a recou- 
pling scheme based on the symmetrized embedded proce- 
dure discussed in Sec. Ill Bl while choosing r according to 
the modified condition ((14|) . we notice that the resulting 
fidelity (squares) is fitted well by a function of the form 
cxp(— Crs/n-swhh) with c rs ~ 0.41. If symmetrization is 
omitted an intermediate behavior is obtained (circles). 



IV. NUMERICAL SIMULATION OF A 
QUANTUM ALGORITHM 

In this section the question is explored how much can 
be gained by stabilizing an iterative quantum algorithm 
by the embedded recoupling scheme of Sec. Ill Bl 



A. Quantum computation with the gate 

For purposes of quantum computation one needs to 
know how to perform two-qubit entanglement gates, such 
as the controlled-not-gate (CNot-gate) or the controlled- 
phase-gate (CP (y;)-gate), on the basis of the recoupled 
Hamiltonian H 1 ^ . Definitely, such quantum gates can be 
performed only between qubits k and I which are coupled, 
i. e. for which Jf-i ^ 0. Therfore, in order to be able to 
entangle any two qubits of a quantum computer it is nec- 
essary to swap qubit pairs with vanishing coupling con- 
stants to neighboring positions. Fortunately, such a uni- 
tary swapping gate can be realized easily by the unitary 
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FIG. 4: Quantum gates of the Swap, CNot, and controlled-phase gate (CP (ip)): R± a (<fi) = exp(=pia Q1 /2 • ip), a £ {x,y,z}. 
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FIG. 5: Left: The qubits of the nine-qubit quantum informa- 
tion processor are arranged on a lattice. The lines connecting 
qubits i and j indicate the values of the coupling constants 
Jij\ right: The table shows the two different values of the 

(2) 

coupling constants Jjy for each qubit pair (k,l). 



phase gate of Eq. (JHJ) because Swap fc; 



U tv/4- 



Through- 



out the rest of this paper we will use the quantum phase 
gate U™ a as a basic building block for all two-qubit quan- 
tum gates. Thus, the quantum Swap fcr gate consists of 
the repeated application of two such gates. For the real- 
ization of other two-qubit quantum gates repeated appli- 
cations of this /g-gate in combination with single-qubit 
gates are required. In Fig. 0] basic gate decompositions 
are depicted for the CNot-gate, the CP (ip), and for the 
Swap-gate. These decompositions will be used in the next 
section for the simulation of a quantum algorithm. The 
^/g-gate itself can be generated approximately either by 
repeated application of the original or of the randomized 
Super- WHH recoupling sequence using either condition 
© or relation (|14|t for the determination of the free evo- 
lution time t between successive fast pulses. 



equal while the coupling constants of diagonal neighbors 
are smaller by a factor of 2~ 3 / 2 due to the larger distance 
between them. Non-neighboring qubits are assumed to 
be uncoupled. According to relation l|14|l this implies 
that in the embedded recoupling scheme two different 
time intervals t are required for the free evolutions. The 

(2) 

values of the coupling strengths for the 9-qubit lat- 
tice used in our subsequent simulation are apparent from 
the table of Fig. 

In the following it is assumed that a quantum algo- 
rithm is performed on this quantum information proces- 
sor according to the following rules: 

1. Single-qubit gates are performed instantaneously 
and perfectly. 

2. Two-qubit gates between vertical or horizontal 
neighboring qubits are performed by repeated ap- 
plications of the unitary flwg-gate of Sec. IIV Al in 
combination with single-qubit gates as illustrated 
in Fig. QJ The fT^/g-gate itself is generated by 
applying Super- WHH sequences n sw hh times as in- 
dicated in Fig. |21 

3. If the target qubits of a two-qubit gate are not ver- 
tical or horizontal neighbors they are moved into 
such positions by applying a sequence of Swap-gates 
according to the following simple strategy [Tg: If 
the vertical position of the qubits is the same, move 
the lower qubit to the upper one. Otherwise, move 
the lower one to the same horizontal position and 
afterwards move the left one as far as necessary to 
the right. 



B. Lattice model of a quantum computer 

For the subsequent numerical simulations of a quan- 
tum algorithm we consider a quantum information pro- 
cessor consisting of n q = 9 qubits which are arranged on 
a lattice as indicated in Fig. [3] The coupling constants 
of vertical or horizontal qubit pairs are assumed to be 



A Super- WHH sequence is always applied in such 
a way that the qubit whose physical position has 
the smaller label (compare with Fig. \5§ is qubit 
k in W£y, i.e. it is transformed by the a x trans- 
formations. The gate sequence of the CP (</?)-gate 
(compare with Fig. HJ is applied in such a way that 
the first single-qubit gate is applied always to the 
qubit with the smaller label. 
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C. The Quantum Algorithm 

In order to investigate the stabilizing properties of the 
embedded recoupling scheme the quantum algorithm of 
the quantum sawtooth map 13] is simulated according 
to the rules of Sec. IIVBI One iteration of the quantum 
sawtooth map transforms an initial n g -qubit quantum 
state |\&) to the quantum state 

\iS,y = e -iTpV2 e -ikv0)/2^ (16) 

with the sawtooth potential V{6) = (9-tt) 2 (0 < 9 < 2tt) 
and the (dimcnsionlcss) momentum operator p. Initially 
the nine-qubit quantum information processor is pre- 
pared in the momentum eigenstate \^>) = |100110011). 
The (dimensionless) parameters of the sawtooth map 
are assumed to have the same values as in the previ- 
ous simulations of Ref. HH, i.e. T = 271-/2™* and 
K = kT = —0.5. Therefore, in Husimi-distributions, 
such as the ones presented in Fig. 03 the dynamics of 
the sawtooth map is restricted to a phase-space cell of 
size 2tt x 2ir and its corresponding classical dynamics are 
integrable. In these Husimi-distributions the initial state 
corresponds to a horizontal line slightly above the middle. 

Our gate decomposition of the quantum algorithm of 
this sawtooth map consists of n g = 2n 2 + 2n q quantum 
gates. It differs slightly from the gate decomposition of 
Ref. [Rj] since we are using the CP (Vhgate instead of 
the four-phases phase gate of Ref. 13]. In particular, 
[2 ■ n q (n q + l)/2] quantum gates originate from the two 
quantum Fourier transforms after which the inversion of 
the qubit positions is taken care of by relabeling instead 
of swapping. 



D. Numerical Results 

In Figs. 16171 and[S]results of our numerical simulations 
of the fidelity 

f(t) = \(V(t)\y id Ut))\ 2 (17) 

are presented for different numbers of repetitions n sw hh £ 
{5,6, ...,10, 16} of the Super- WrlH sequences. For each 
value of n sw hh we calculated the fidelity of the quan- 
tum state of the quantum sawtooth map for up to 
t = 300 iterations as well as the corresponding Husimi- 
distributions. 

The quadratic-in-time fidelity decay of the non- 
randomized recoupling scheme is clearly apparent from 
Figs. El and [7| (The corresponding fidelities are plotted 
in black color). This decay is caused by the coherent ac- 
cumulation of errors due to the second-order AHT-term 
of the Super- WHH sequences involved in the realization 
of the unitary L/^/g-gate ^3- The i-dependence of the 
fidelity can be fitted by the function 

/(t)-exp(-c-i 2 Kwhh) (18) 
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FIG. 6: Upper plot: Fidelity plots of the quantum saw- 
to oth map implemented with the original recoupling scheme 
of lYamaguchi et al.| l f (black) and the corresponding plots of 
the symmetrized embedded recoupling scheme (red): Dashed 
curves show the fidelity estimations according to Eqs. I|18|l 
and 11911 . Lower plot: Fidelity plots of the restricted random- 
ized (blue) and of the symmetrized (red) embedded recoupling 
scheme. 
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FIG. 7: Logarithmic fidelity plots of the quantum sawtooth 
map with n aw hh = 5 (left) and n sw hh = 10 (right): the original 
\Yamaauchi et ql.|l| scheme (black) , the restricted but unsym- 
metrized embedded approach (blue), and the restricted and 
symmetrized embedded scheme (red). 



with c w 0.87 (compare with the seven lowest dashed 
lines of the upper picture of Fig. |HJ). According to Ref. 

there should also be a linear contribution in the ex- 
ponent of (|18|l which dominates the fidelity decay for 
small numbers of iterations. Neglecting this linear con- 
tribution is the reason for the slightly imperfect overlap 
of our fitted fidelities with the corresponding numerical 
results. 

Simulations based on the restricted embedded recou- 
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flswhh = 6 n sw hh = 8 W sw hh = 10 



FIG. 8: Husimi-distributions of the quantum states resulting 
from the qua ntum sawtooth map : (Upper row) The origi- 
nal scheme of lYamaguchi et al.| l f . (Middle row) the embed- 
ded but unsymmetrized approach, (Lower row) the embedded 
and symmetrized scheme. The C/^/g-gates used in the com- 
putations consist of n sw hh = {6, 8, 10} Super- WHH sequences 
(from left to right). These distributions are averaged over 
290 < t < 299 numbers t of iterations of the sawtooth map. 

pling scheme without symmetrization are shown in Fig. 
El (lower part, blue plots) and Fig. (blue plots). The 
fidelity decay is suppressed significantly but on the time 
scale of these plots it is still quadratic in time. This orig- 
inates from the fact that terms of the Hamiltonian of Eq. 
© of the form a^a l a , a £ {x,y, z}, are not eliminated 
by the PAREC method without symmetrization. 

Using the randomized and symmetrized Super- WHH 
sequence together with the appropriately chosen free evo- 
lution times (compare with Eq. I|14|l ) it is possible to get 
an almost linear-in-time fidelity decay at least on time 
scales where errors of the order of 0(J( Jr/fi) 4 ) are neg- 
ligible (compare with Figs. El and [7| (red plots)). In these 
cases the fidelity decay can be fitted by the function 

f(t) = exp(-c rs • t/n 5 swhh ) (19) 

with c rs w 7.85 (compare with the six upper dashed lines 
of the upper picture of Fig. which are almost indistin- 
guishable from the corresponding full curves). 

V. CONCLUSIONS 

We showed how a selective recoupling scheme can be 
embedded into a stochastic decoupling scheme in such a 



way that a particularly selected coupling is achieved and 
that, in addition, the coherent accumulation of higher- 
order errors is suppressed significantly. Even if compu- 
tation times of a quantum information processor are so 
long that the residual higher-order interaction term of 
Eq. Ijlijfl of the order or 0( J( Jr/fi.) 4 ) is no longer negligi- 
ble, it is possible to suppress also these errors significantly 
by a suitable adjustment of the free evolution time r in- 
volved in the realization of the relevant two-qubit gates 
(t^,-gates). In generalization of the procedure discussed 
in Sec. Ill Bl (compare with Eq. I|14l) ) this can be achieved 
either by explicitly calculating the fourth-order contribu- 
tion of AHT and by solving the corresponding implicit 
equation of fifth order for r involving renormalized cou- 
pling strengths or, alternatively, adjusting the value of r 
so that the resulting fidelity decay is as small as possible. 

Basic properties of our embedded stabilization scheme 
were analyzed for a single two-qubit gate. In particular, 
it was demonstrated that our proposed embedded sym- 
metrized recoupling scheme results in an improvement of 
the scaling of the error of a swapping gate with n~^ hh 
instead of ti~ j^, ■ Thereby, n sw hh denotes the number of 
repetitions of an embedded Super- WHH sequence which 
are required for the realization of the phase gate. There- 
fore, in our embedded recoupling scheme fewer numbers 
of repetitions of Super- WHH sequences are necessary for 
achieving a particular degree of error suppression. Typ- 
ically, this also implies fewer pulses which are required 
for performing a quantum computation with a particular 
error tolerance. This aspect is apparent from the upper 
plot of Fig. Elwhere at t ~ 70 iterations the fidelity of the 
original recoupling scheme with n sw hh — 16 is the same as 
the one of the embedded symmetrized recoupling scheme 
with n sw hh = 6. 

Finally, we want to address briefly the additional phase 
shifts each selective pulse is accompanied by in any re- 
coupling scheme. These phase shifts (compare, e. g., with 
Eq. (24) of Ref. 0) originate from the transformation 
from the interaction picture to the Schrodinger picture. If 
necessary, such a phase shift can be eliminated by letting 
the system evolve freely after the application of a selec- 
tive pulse for the same amount of time thereby embed- 
ding this additional time evolution between two broad- 
band 7T pulses. This way an additional phase shift with 
the opposite sign is generated. 
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